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A STUDT 0? THE DIEFUSE GALACTIC GAIMA RADIATION 


Carl Tlchtel and Donald Knlffen 
NASA/Goddard Space Flight Center 
Greenbelt, MD 20771 


I. INTRODUCTION 

Assuming cosmic rays pervade the Galaxy, they necessarily produced high energy 
T-rays as they Interact with the interstellar matter and photons. The cosmic 
ray nucleon Interactions give rise to y rays primarily through the decay cf x" 
mesons, giving a unique spectrum with a maximum at approximately 6R MeV. 

Cosmic ray electrons produce y rays through hremsstrahlung, hut with a 
markedly different energy spectral shape, one which decreases monotonlcally 
with energy. Cosmic ray electrons also Interact with the Interstellar 
starlight, optical and infrared photons, and the blackbody radiation through 
the Compton process. Finally, cosmic ray electrons can Interact with magnetic 
fields giving rise to synchrotron or curvature radiation, but these processes 
are much less Important then the others previously mentioned for the galactic 
diffuse radiation and will not be discussed here. 

Extensive work has already been performed on the calculation of the source 
functions for these various y radiations and the intensity to be expected in 
the vicinity of the eular system. (For a general review see Chapter 5 of 
Flchtel and Trombka, 1981.) However, several recent devlopments make a 
reexamination and extension of this work worthwhile. These Include the 
recently published detailed results of high energy galactic y-radlatlon 
obtained with the COS-B satellite (Mayer-Hasselwander et al. , 1982), further 
evaluations of the 21 cm absorption radiation In the galaxy and, hence, the 
atomic hydrogen density In the Inner galaxy, considerations related to the 
molecular hydrogen density normalisation, recent results on molecular clouds, 
and improved theoretical calculations on the nucleon-nucleon source function. 
The developments related to the galactic matter distribution and y-ray 
production will be considered In the next section and Incorporated Into the 
general Y~ray production calculation. The predictions will then be compared 
to the recently published COS-B high energy Y“fay results in section III. 

TI. DIFFUSE GALACTIC GAMMA RAY PRODUCTION 

(a) Galactic Matter, Photon, and Cosmic Ray Distributions 

With regard to the matter, the relevant concern is the galactic diffuse matter 
In the form of atoms, molecules, ions and dust. The latter two are believed 
to be minor constituents and, hence, unimportant for Y~ray production through 
cosmic ray interactions. Hydrogen Is the primary component of both the atomic 
and molecular matter. Helium and heavy nuclei add about 55X more to the y~ray 
production. It is assumed these latter nuclei huve a distribution In the 
galaxy similar to hydrogen, although little Is known about them. Both atomic 
and molecular hyd*'ogen are known to be confined to a narrow disk 0.12 kps 
in scale height for atomic hydrogen, but see specifically Knlffen, Flchtel and 
Thompson, 1977, for the model of scale heights used here) with the molecular 
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hydrcgen distribution apparently somewhat narrower (e.{ , Gordon and Burton, 
1976; Solomon and Sanders, 1980). Atomic hyt’rogen revt Is Its presence 
through the emission of the 21 cm line; however, there smalns some 
uncertainty in the density of the inner galactic region related to the 
absorption correction. Recent work (e.g* Dickey et al., 1982; Thaddeus, 1982) 
suggests that the absorption had previously been somewhat underestimated and 
that the density In the region of 4 to 5 kps from the galactic center may be 
greater than previously estimated perhaps by a factor of In this work the 

atomic hydrogen density distribution as a function of radius from the galactic 
center of Gordon and Burton (1976) was used, but modified so that the atosdc 
hydrogen density In the Inner region was Increased by a factor of 1.5, and the 
closer densities were Increased less In accordance with the amount of 
Intervening matter. This density Is also modulated for the galactic arms In a 
manner to be described later. 

The density distribution of molecular hydrogen Is measured less directly. At 
present, the best approach appears to continue to be through the observations 
of the 2.6 mm spectral line “)f which the distribution of cold 

Interstellar matter is Inferred. The nature of the Interpretation of these 
measurements In terms of what Is really desired makes the otolecular hydrogen 
density distribution less certain than that of the atomic hydrogen. The 
average galactic radial distribution of molecular and atomic hydrogen deduced 
by Gordon and Burton (1976) does show reasonably clearly that the molecular 
hydrogen to atomic hydrogen ratio Is larger In the Inner galaxy than It Is In 
the outer galaxy even If the absolute Intenlsty of molecular hydrogen Is still 
quite uncertain. For the work here, the molecular hydrogen density 
normalization Is treated as an adjustlble parameter In the range from that 
estimated by Gordon and Burton (1976) to a factor of 0.4 smaller. The final 
value actually used was 0.6 smaller. It should be noted that the CO 
observations Indicated that the great majority of the molecular hydrogen is In 
clouds. The recent work of Solomon and Sanders (1980) has. In fact, auggested 
that the Interstellar medium Is dominated by massive cloud complexes. 

Although the translation of the observations into a galactic spaclal 
distribution Is difficult, on a broad scale the density profile Is reasonably 
well known, even though details of arm structure are not always agreed on by 
all workers In the field. A general spiral pattern does appear to emerge. In 
addition to the 21 cm data the distribution of continuum radiation (I-andecker 
and Wleleblnski, 1970; Price, 1974), y radiation (Blgnaml, et al., 1975), Hill 
regions (Georgelln and Georgelln, 1976), supernova remnants (Clark and 
Caswell, 1976), pulsars (Selradakls, 1976), and Infrared emission (Hayakawa et 
al., 1976) are all consistent with the existence of spiral structure In the 
galaxy. Until recently, it had not been clear whether molecular clouds were 
associated with spiral structure. However, now on the basis of a high sample 
survey and observations in both the first and second quadrants of the galactic 
plane, Cohen et al. (1980) have shown the existence of the molecular 
counterparts of the five classical 21 cm spiral rams segments In these 
quadrants, namely the Perseus arm, the Local arm, the Satlgarlus arm, the 
Scutum arm, and the 4 kpc arm. Kutner and Mead ^981) have even Identified 
arms through CO measurements In the outer galaxy. The specific spiral pattern 
that will be used here Is that of Georgelln and Georgelln (1976), and the 
matter density In the plane will be assumed to follow the pattern described by 
Knlffen and Flchtel (1981) relative to the am to Interam density ratio and 
the am width. 
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ORIGJMU PACI f.3 

OF POOR QUALITY 


For the photon dlstributlone Kniffen end Flehtel (1981) neing reenlte of Boeel 
^t el. (1982) on the Infrered voluae eaieelvlty end e aodel of Beheell end 
Sonelre (1980) for the eterll :ht distribution, obtelned photon densities end, 
hence, source function for the Conpton salsslon es e function of position In 
the gelexy. These will be used here. 

With regerd to the cosalc rey distribution in the gelexy (see pertleulerly 
Kniffen end Flehtel, 1981, end Flehtel et el., 1976), It will be essuaed thet 
the nucleonic cosalc rey composition end energy spectrum remeln unehenged 
throughout the gelexy end thet the electron spectnn chenges only In e second 
order menner as the density chenges. Again, for the reasons described In the 
above works, the cosalc rey density In the plane will be essuaed to be 
proportional to the matter density on the scale of arms and, perpendicular to 
the plane, to have a gausslan distribution with a scale high of 0.6 kpc. This 
Irtter number Is based on the radio continuum measurements of Cane (1977) and 
the assumption that the galactic magnetic fields energy density and the cosmic 
ray energy density have the same scale height. This scale height for the 
cosmic rays Is somewhat less than that used previously, and the primary effect 
is a relative reduction In the Compton contribution. 

(b) Gamma Ray Source Function and Calculation of Prediced Intensities 


The detailed calculations asscociated with the pro action of energetic y rays 
through cosmic ray nucleons interaction with Interstellar matter including all 
the primary cosmic ray and Interstellar matter components, all the secondaries 
and their decay products, the angular distribution, and the energy spectrum 
are very detailed and lengthly. These calculations have, however, been 
performed. See, for example, Cavallo and Gould (1971), Stacker (1971), Badwar 
and Stephens (1977), and Morris (1982). The latter's work Is based on the 
inclusion of substantial recent experimental work Into a model of y~ray 
production which rerulres momenttim balance in the center of mass system. The 
predicted spectra are shown in Figure 1. The significance of the differences 
will become more apparent in the discussion of the experimental results in the 
next section. 


Fig. 1: Spectrum of the 

gamma-ray source function for 
cosmic ray nucleon, 
interstellar nucleon 
interatlons as a function of 
energy (Mrrrls, 1982). 
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OF POOR QUA-UY 

The cosmic ray electron* matter Y^tay production can be calculated ualng the 
bramsstrahlung croaa-aeetlon formulas of Koch and Kota (1959). The predicted 
radiation in the region below about 10^ MeV la uncertain even locally In our 
galaxy because the Interstellar cosmic ray electron spectrum Is not wall known 
at low energies where the electron spectrum observed near the Earth has 
undergone strong solar modulation. For the work here* the calculated source 
functions of Flchtel et al. (1976) based on the cross-sections of Koch and 
Mots (1959) will be used here. 

The calculations associated with the production of Compton Y rays have been 
performed in some detail for the cases cf astrophysical Interest by Ctnsburg 
and Syrovatskll (1965). Cosmic ray electrons interact with galactic starlight 
photons* for which the optical and Infrared ranges are the Important ones* and 
with the universal blackbody radiation. The source functions of these 
Interactions are much smaller In the galactic plane than that for 
bremsstrahlung. However, the total contribution to the galactic y radiation 
Is significant because the cosmic ray and stellar photon scale height above 
the galactic plane are gre«iter than those of the matter. The results 
presented here for the Compton radiation will be based on the source functions 
of Knlffen and Flchtel (1981). 

The Intensities In any direction are then calculated In a manner described, 
for example* by Flchtel and Trombka (1981) with consideration of the angular 
resolution of the Instrument being taking Into account where appropriate. The 
predicted variation with longlt\«de for the various Compton components and the 
cosmic ray* matter Interaction component are shown In Figure 2. 


Fig. 2: Gamma ray photon 

Intensity as a function of 
galactic latitude predicted by 
the model discussed here for 
the Compton components and the 
total cosmic ray* matter 
Interaction component. 
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?:il. GAMMA RAY RESULTS AND THEIR INTERPRETATION 

Based on the source functions* the natter dlstrlhutlons* and the assumptions 
just discussed, the expected Y**Tay Intensity has been calculated. The 
prediction of the T~ray Intensity vlth the molecular hydrogen assined to be 
0.6 that of Gordon and Burton (1976) Is compared to the SAS-2 and COS-B 
longitude distributions In Figures 3 and 4 and to the energy spectrum In the 
galactic center region in Figure 5. Considering the uncertainty In the point 
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Fig. 3: The high energy E > 100 MeV y-ray intensity as a 

function of longitude for -10* < b < 10* from the SAS-2 data 
(Hartman et al., 1979) compared to the model discussed here. 


source contribution aad the mass distribution, the agreement between the data 
and the predlcod curves seems reasonably good. Regarding Figures 4 and 5, 
there are two comments. First, If the older cosmic ray nucleon source 
function of Stacker (1971) had been used rather than that of Morris (1982), 
the agreement in the 300 MeV to 5000 MeV energy interval would be quite poor 
with the theory predicting 1.7 to 2 times as many y reys in the center than 
reported. If this older source function were correct, then either the cosmic 
ray nucleons play a relatively small n 'e for some reason, or the measured 
high energy COS-B y-ray intensity is too low. Even with the more recent 
nucleon-nucleon calculations, the >300 Me\ intensity seen by COS-n is 
generally a bit lower than expected relati' o to the 70-150 MeV intensity. The 
most straightforward adjustment to the model to account for this poslble 
difficulty would be an enhancement of the electron component coupled with a 
decrease in the normalization parameters for molecular hydrogen. The second 
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SAUaiC LONGITUDE i 

Fig. 4: Gamma ray Intensity as a function of longitude averaged over the 

latitude range -10" < b < 10" from 70 MeV - 150 MeV, 150 MeV - 300 MeV. and 
300 NeV * 5000 MeV from the COS-B data (Mayer-Hasselwander et al., 1982) 
compared to the model discussed here shown by the solid line and a constant 
cosmic ray density distribution model shown by the dashed line. 


Fig. 5: Energy spectrum of 

the galactic Y radiation for a 
region near the galactic 
center. The calculated 
spectra are based on the work 
of Knlffen and Flchtel, (1981) 
modified to Include the recent 
nucleon-nucleon calculations 
of Morris (1982). The COS-B 
data are those of Mayer- 
Hasselwander et al. (1982), 
and the SAS-2 data are those 
of Hartman et al. (1979). 
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cotment is that, although tha Intensitlas for tha cantar and tha aoureas 
(195,5), PSR 0531+21, PSR 0833-45 naasurad with tha SAS-2 and COS-B 
Instrumants ara In good agraanant with each othar, tha COS-B instruacnt 
appaars to saa a somawhat largar (about 25Z) anticantar diffusa intensity. 

The SAS-2 instnimant, which sees the lower diffuse intensity, had essentially 
no background; the vary small extragalactic Isotropic background is not 
substracted in Figure 1. In the case of COS-B, which seas tha higher diffusa 
intensity, the background estimated by tha COS-B collaboration has bean 
subtracted. 

Tiic latitude distribution predicted by the model taking into account the 
angular resolution of the COS-B instrument generally agrees well with the COS- 
B observations up to latitudes of about eight degrees. Beyond this point the 
data tends to exceed the predictions of the model. An alternative description 
is that there is a small constant difference at all latitudes. This 
difference is consistent with the apparent excess in the COS-B dati^ deduced 
from the longitude distribution just discussed. 

A constant cosmic ray density, as might be predicted in the universal cosmic 
ray model, predicts a rather small T^ray Intensity from the central region. 

The dashed curve in Figure 4 refers to a model identical to the one juct 
discussed except that the cosmic rays are assumed to be constant. Hence, if 
contrary to theoretical expectations, the cosmic ray density should be 
uniform, a relativel> large point source contribution %iould be needed in the 
galactic center region. 

It should be mentioned that there is also an unresolved point source 
contribution to the "diffuse" radiation measured by the SAS-2 and COS-B y-ray 
instruments since the limited angular resolution of these instruments does not 
permit the separation of point sources. It is quite difficult to estimate 
this contribution; however, several factors suggest that point sources may . ot 
be a major contributor (see, for example, Cesarsky, 1980). These include the 
uniformity of the energy spectrum just discussed and, as will be seen, the y- 
ray luminosity of the galaxy and its distribution being about what would be 
expected from the diffuse sources. For the purpose of this paper, the reader 
is simply asked to keep in mind that there is some point source contribution 
yet to be determined which at least for the moment is assxmed to be small, but 
not zero. 

IV. CONCLUDING REMARKS 

Considering the uncertainty in the point source contributions and the 
differences between the SAS-2 and COS-B data due presumably to the more 
limited statistics of the SAS-2 data and the need to subtract a substantial 
and necessarily somewhat uncertain background from the COS-B data, the 
agreement between the theoretical predictions and the y-ray data seems quite 
reasonable. It should be noted that, in general, considering the difficulties 
and pioneering nature of the experiments, the agreement between the SAS-2 and 
COS-B data is remarkably good In terms of general intensity level, energy 
spectra, and relative distribution. At present, detailed refinements to the 
model being discussed here seem inappropriate: however, when more detailed 
information exists, as it may soon, on the location of the large molecular 
clouds and especially the close ones, a further pursuit of the problem might 
be worthwhile. 
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